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Abstract: In 2004, Creminelli and Zaldarriaga proposed a consistency relation for the pri- 
mordial curvature perturbation of all single-field inflation models; it related the bispectrum 
in the squeezed limit to the spectral tilt. We have developed a technique, based in part on 
the Creminelli and Zaldarriaga argument, that can greatly simplify the calculation of the 
squeezed-limit bispectrum using the in-in formalism; we were able to arrive at a generic 
formula that does not rely on a slow-roll approximation. Using our formula, we explicitly 
tested the consistency relation for power-law inflation and for an exactly scale-invariant 
model by Starobinsky; for the latter model, Creminelli and Zaldarriaga's argument pre- 
dicts a vanishing bispectrum whereas our quantum calculation shows a non-zero bispectrum 
that approaches zero in the long-wavelength limit and for inflation with a large number of 
e-folds. 
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1. Introduction 

Inflation was originally proposed in the 1980's to solve the monopole problem, the horizon 
problem, and the flatness problem [Q, |2|, ^, |^, ^. However, its strongest validation has come 
from the study of cosmic perturbations, since inflation provides a natural mechanism for 
producing the large scale perturbations observed in the CMB and in large-scale structure 
1^, ^, 10, 11 1 . We now know that the fluctuation produced were nearly scale invariant 



and were of the order of 1 part in 10^ of the energy density. 

The observable consequences of inflation are usually written on terms of correlation 
functions. One of the most promising avenues for further discrimination among models 
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of inflation is the study of non-Gaussianity - that is, the connected part of the three- 
point or higher-point correlation functions - in the primordial curvature perturbation ^. 
C is conveniently defined as the scalar perturbation in comoving gauge^ , where the metric 
including only scalar perturbations is 

ds^ = -(1 + A)dt^ + 2aBidx' dt + {t)e^^^''Ux'^ . (1.1) 

At later stages, Q (or nearly equivalently, the Newtonian potential $) can be used to 
predict perturbations in large-scale structure or the cosmic microwave background. For 
single-field inflation, Q is particularly useful because it is conserved outside the horizon 
|1C, 13, 14]. The lowest-order non-Gaussianity is the three-point function or bispectrum, 
which we parametrize as follows: 

(CkiCk^Cka) = (27r)3<5{3)(^k,) i?^(ki,k2,k3). (1.2) 

For future reference, we will similarly write the power spectrum as 

(Ck.Ck^) = (2vr)=^5=^ (ki + k2) Pi^{k), (1.3) 

and define the spectral tilt as — 1 = so that we have P oc A;"'^"'""'^. 

The first results about primordial non-Gaussianity were in terms of the non-linearity 
parameter /nl (sometimes called /^l); /nl remains the best-constrained non-Gaussianity 
measurement, with /nl = 32 zt 21. /nl parametrizes the part of the bispectrum that has 
the form 

(CkiCk^Cka) = {2TTf5^%Y>r) pNL[Pdh)Pdk2)+Pdk2)Pdh)+Pdk3)Pdkl)] . 

Since we measure that Pdk) oc k^"^ [|l^], it is easy to see that this bispectrum has a squeezed 
"shaped" — i.e. that it peeks in the squeezed limit — and we can calculate its squeezed 
limit form^ 

(Ck,Ck2Ck3)fc3«,,,, = (2vr)3 5(3)(^k,) ^f^^P^(k^)Pdh). (1.4) 

Clearly, an important observational test for any potential model of inflation is that its 
predicted /nl be compatible with observed results. 

In 2004, Creminelli and Zaldarriaga |17|, generalizing an observation from an earlier 
paper by Maldacena [|l4|, used a clever argument to impose a consistency relation on the 
three-point function of single-field inflation, or more specifically, on the three-point function 
in the local limit (i.e. when one of the wavenumbers is considerably smaller than the other 
two). They demonstrated, using a purely classical argument, that for any single- field 
inflation model: 

(Ck.Ck2Ck3)fc3«;i„fe = (2vr)=^5(3)(^ki) {l-ns)Pdh)Pdk3). (1-5) 



^This corresponds to a gauge where Ht ~ (in the notation of [^2|) and, for single-field inflation, 
5(p = 0. 



Remember that, in the squeezed limit fca <C fci so that ki ~ h 
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In other words, for any single-field model, the consistency relation tells us 



/NL = ^(l-n,). (1.6) 

Since 1 — ns = 0.037 ||l6[], every single-field model should produce /nl ~ 0.02; it therefore 
follows that, according to the consistency relation, any measurement of /nl of order 1 or 
higher rules out single-field inflation. Thus, the consistency relation can potentially help 
rule out a very general class of models. For this reason, it is important to thoroughly 
understand it. 

In particular, it is interesting to consider what can be learned from trying to verify the 
consistency relation using a quantum field theory approach rather than the classical meth- 
ods of the Creminelli and Zaldarriaga proof. One can see one fruit of this approach when 
we test the Starobinsky model in subsection where we can quantitatively determine 
corrections to the consistency relation prediction in the case of a finite inflationary period 
where all modes leave the horizon. 

In this paper, we will demonstrate a technique for calculating the bispectrum in the 
local limit using the in-in formalism without first calculating the general bispectrum; we 
will then verify the consistency relation for three models. In section we will review the 
argument in [^], since it shares some assumptions with our technique. In section we will 
outline our approach. In the next two sections, we will develop our approach in more detail. 
First, in section ^, we derive the action we will use to calculate the bispectrum. Then, in 
section ^ we perform steps outlined earlier and arrive at ( 5.1^ ), the general formula for the 



bispectrum. In section y, we use our formalism to explicitly verify the consistency relation 
for slow-roll inflation (a known result), power-law inflation and an exactly solvable model 
jl8U by Starobinsky, the latter two with no slow-roll approximation. Finally, in section ^ 
we will discuss the implications of our results as well as mention other situations where 
our formalism would be useful. In appendix we outline some of the manipulations in 
section || more explicitly; in appendix ^ we list the properties of Hankel functions used in 



subsection 6.2 on power-law inflation. 



2. Review of the Creminelli and Zaldarriaga argument 

We are interested in calculating (CkiCk2Ck3) in the limit that k-^ <^ ki, k2 (so that ki ~ k2). 
Mathematically, 



(Ck,Ck,Ck3) = ((Ck,Ck,>^^^Ck3), (2.1) 

where we deflne ^ . . . V to be the expectation value of . . . given that Cka has a partic- 
ular value. Our focus (both here and in later sections) will revolve around calculating 
(CkiCka)^^^- 

We will evaluate (CkiCk^)/- after the ki,k2 modes have crossed the horizon so that 
the ks mode will have crossed the horizon in the distant past. Thus, will be part of 
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an essentially classical [|T9| , 20 j background which affects the scalar field through the 
metric. Considering only modes far outside the horizon, the metric is 



(2.2) 



where 



C^(x) 



(Note that this also fixes our Fourier convention). 

For the next step, let us consider the equivalent real space correlation ^C^)^^ (xi, X2) 
(we will use ( . . instead of ( . . . )^ for the remainder of this section^). We know 

that the background perturbation (^^ is small, so it makes sense to expand the correlation 
function in a power series about and keep only the first term: 

(C^>^,(x,Ax) = (C'>o(Ax) + I d'kC^ (^|^^_^(C2>^,(x,Ax)^ +... , (2.3) 

where x = (xi + X2) /2 and Ax = X2 — xi. 

To evaluate this expression, we need to know more about (^C^^^g. For ks small enough, 

we can perform a coordinate change x — t- x' = e^^^^^x to put ( p.2[ ) in the form of the 
unperturbed FRW metric. In these new coordinates, the background is unperturbed so 



<C^>.,(x,Ax) « (C^>o(|x'2 -x;|) « (C^>^(e^"WAx). 



(2.4) 



Thus, 



_S_ 



(C^) (xAx)-^lMo^M 



(2.5) 



substituting this into (^), moving to Fourier space, and correlating with Cka^ 



(CkiCkaCka) = \ (CkiCk2)(^^^Ck3 



din Ax 



(%) 



.-(..)3..3.(E.,)m)^[fM 

= (27r)3 5(3)(^k,)P(A;i)P(A;3)(l-n5), 
thus yielding the result of 0. Note that ks = (ki - k2)/2 ki. 



Otherwise, (2.3) becomes very confusing. 



-4- 



3. Outline of our formalism 



In this section, we outline our approach, which we then follow in more detail over the next 2 
sections. As in the proof in |17|, we will calculate (C^)^|^ ^^'^ then use = (('^^)ck 



<3l 



We will also split ( into 2 parts: a large-scale, classical, background part Cl and a small- 
scale part which undergoes quantum fluctuations; we let 



Cl - / l^C^e'^-^ Cs = TTT-^Cke^'^-, (3.1) 



(where we choose k^: such that < k^: <^ ki , , so that = Cl + Cs- Note that our method 



is distinct from stochastic inflation |21| because we use canonical quantization. Also, we do 
not worry about the precise value of the cutoff fc^, (an issue which takes on some significance 
in Stochastic inflation; see, for example p^) because we are performing our calculation 



near the time that ki,k2 exit the horizon, so that Cl (which has k < k^ <^ ki) is almost 



completely classical pO[ . 

To calculate (C^)^j^ accurately, we will use cosmological perturbation theory. We start 
with the Lagrangian for a single scalar-field (p and we work in a flat Friedman-Robertson- 
Walker (FRW) background. By an appropriate gauge choice, we can take C to be our degree 
of freedom (instead of (p); C is more convenient for calculations since it is conserved outside 
the horizon. Since we ultimately want the three-point vertex {C^), we need the Lagrangian 
up to third order in Then, we split ( in the Lagrangian, using C = Cl + Cs, and keep 
terms of order C| and C|Cl- Also, we will only need the terms that are lowest order in 
spacial derivatives and in time derivatives, since these will dominate in the squeezed limit; 
this will allow us to drop a number of terms. If the space and time derivatives on a term 
act in opposite ways it may not be clear how such a term compares to other terms and we 
have to keep it^. 

We can now canonically quantize Cs- We treat the terms of order Cs as providing 
the equation of motion for the mode functions of Cs and we consider the terms of order 
CsCl as being perturbations. Then, we use the in-in formalism (a technique in quantum 
perturbation theory useful in cosmology) to calculate (Ci)^j^ • Finally, we correlate with 

Cka and arrive at (C^)- 

From this point forward, we will sometimes omit the subscript on (. . .)^^ , except where 
it adds clarity. 

4. Deriving the action 



In this section, we will essentially follow |14|. 



We will work in units where c = h = 1 and rup'^ = SttGn = 1- We assume that, on 
average, inflation takes place on a homogeneous, nearly flat background; thus, we have the 



*For example, if we compare a term like d~^(, (where is the inverse Laplacian) with the first 
term has 1 time derivative but -2 space derivatives while the second has time derivatives and space 
derivatives. Thus, we can not say which one dominates and we need to keep both. This problem actually 
only comes up once and it turns out we can safely ignore that term. 
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following FRW background metric: 



ds^ = -dt-^ + a\t)dx\ (4.1) 

It is sometimes convenient to use conformal time r/, where r] = J dta~^ so that dt = adrj. 
Dots over characters (e.g (p) refer to derivatives with respect to physical time t, while primes 
(e.g (p') refer to derivatives with respect to conformal time. We define the usual Hubble 
parameter as H = d/a. 

We suppose that a single scalar- field (p with a canonical kinetic term is the inflaton^. 
The action for (p is 

S = ll d^^gx[R - (V<pf - 2V{<P)]. 

We write (/'(x) = (pQ{t) + (5(/)(x), where (pQ is the homogeneous background part of the field 
and 6(p is a small perturbation. Then, from equations of motion and conservation of energy, 
the background (pQ obeys the following: 

2,H^ = lcPo^ + V{<Po), (4.2a) 
H = -]^cPo\ (4.2b) 
O = 0o + 3i?0o + ^'(</>o). (4.2c) 

To determine the dynamics of perturbations, we need to calculate the action for 6(p (or 
equivalently, for ("). The most efficient way of doing this is to use the Arnowitt-Deser-Misner 
formalism |24|, where we write the metric as 

ds'^ = -N'^dt^ + hij{dx' + N'dt){dx^ + N^dt). 

The purpose of the ADM formalism is that we can treat the lapse, shift variables N,N^, 
respectively, as Lagrange multipliers in order to integrate out their degrees of freedom. 
Using these variables, we can rewrite the action as 

S=\ j d^xVh NR^^^ - 2NV + N'^iEijE'^ - E^) + N-^{<i) - N'diCpf - NW^di(pdj(p 
where 

Eij = ^{hij -V,Nj -VjNi), 
E = El 

where the covariant derivatives on Ni are taken with respect to the three-geometry. 

We will adopt comoving gauge 13 1 where, considering only scalar perturbations, 



^The generalization to a non-standard kinetic term is done in 
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In this gauge, we have shifted the degree of freedom from 54> to C is the same quantity 
introduced early, namely the Bardeen variable that is conserved outside the horizon |2^, 26 1. 

We can now exploit the power of the ADM formalism. We find the equations of motion 
for and N, solve them to the necessary order in (2'^'^ order for our purposes), and 
then replace them in the action^. 

Up to third order in we find 



S : 

S2-- 
S3 



where 



S2 + Ss, 



1 i4 



1 00 14>1 

2 </,oF 4 m 



-2 

4 IP 



1 1 

2li 



d-^didj{dixdjC) 



4 i/2 
6L 



(4.3) 



6L 



SS2 



to lowest order, the equation of motion for is 



0. 



This is the action we want. Note that we have not made any slow-roll assumptions. 



5. Calculating the local-limit bispectrum 

As described in section^ our goal is essentially to determine (CkiCiL2)r = (C5kiC5k2)/- ; 
that is, the two-point function calculated with a classical background C^, where = (^kg. 
Ordinarily — i.e. if we do not fix a perturbed background — (CkiCk2) ^(^i + k2) = if 
k2 7^ — ki. Thus, the two-point function {^(^^) is only non-zero insofar as it is affected by 
a contribution from the background, which arises from the third-order term in ([4.3|). We 
will see that the contribution has two pieces: one from the field redefinition of (s in ( ^-31 ) 
and one from the in-in formalism. 

® Again, this is done more explicitly in ]14|. 
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5.1 Finding the action in terms of (^s and (l 

Using C = + C-L) we find the terms of order in 52, which we call So, and the terms 
of order CsCl in 5*3, which we call 5int^- 5*0 will determine the equation of motion for (^5, 
and 5int will be the perturbation to 5*0 and will enter into the in-in formalism to calculate 

From ^2, we see that 



1 



^0 = ;^ / d'^xj^ia-X's-aidCsn 



(5.1) 



We can similarly pick out the C5C-L terms in 53. Note first that the two terms with 
coefficients of have derivatives on all the C's and thus must yield a term with a derivative 
on (^l; these two terms will thus be negligible. From the remaining terms, we substitute in 
for ( and pick out terms with the fewest derivatives on (l^, yielding 



1 



14 



1 



4 16 if6 



aXsdiCsdid-'CL+ 



16iJ6 



'dt 



^0 ^1_ 



2 <PoH 4 



CsCs 



/(C) 



SCs 



(5.2) 



where 

/(C) = 



The terms with a factor of 5Lq/5C,s are best removed by a field redefinition: 



+ 



4>qH 2 m 



ClCn + 



(5.3) 



^There are terms of order CsCi- which could, in principle, be significant. However, after an integration 
by parts, they become equal to 



^4 r 



C=Cl 



Since C^l is a classical field it obeys its equations of motion, so |^ 
contribute. 



= and these terms do not 



C=Cl 



Any terms of order CsCl can also be ignored because they can not contribute to a two-point correlation 

(cl> 

*As mentioned in section H, we retain terms where the effect of the derivatives is unclear, as in the 
did~^C,L term. Actually, it turns out that this term does not contribute, though this is not immediately 
obvious. The reasons are mentioned later in this subsection and in appendix 
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where, in the second hne, 1) we have dropped terms that vanish outside the horizon, since 
the terms in the field redefinition will only be evaluated outside the horizon (as we will 
see in subsection ^.21)^; and 2) we have replaced Cs with Cn since they are effectively equal 
when they multiply terms of order C^. 

Thus, after the redefinition, our action Si^t in terms of Cat becomes 



int 



2H^ 



'^a^CNdiCNdid-''CL+ 



d / 1 



yo _^ 1 



CnCn 



(5.4) 

5.2 Field redefinition 

C,s, rather than (^jsr, is what will remain constant outside the horizon and what we correlate 
with so we need to calculate (C5,kiC5,k2)- To understand the effect of a field redefinition, 
consider a simple case, where Cs = Cn + k and k <^ (n- then 

{Cs,Ut) Cs,Ut)) - (Civ,ki(t) CN,Ut)) + 2(K(t) CM*))- 

We actually have (|5.3D, where Csi^) ~ Cn{x) + f3{t)(L{x)CN{x); in this case, we find"^^ 
(CWt)CWi)> ~ (CWO W0> +2/3(t)CL,ki+k2(0 FT [{C^)^{Ax,t)]{h) 

(this is worked out in appendix A.l). Note that FT [(^^^^(Ax, t)] (fci) = jn^j (t)|^, where 
is the mode function defined in the following section. If we measure at a time t ^ t^,, where 
t^, is the time that ki , /c2 cross the horizon, then FT [(C2)q(Ax, t)] (fc) = \uk{t)\^ = P{k). 
Thus, 

{CsM{t) CsMit)) « {CN,Ut) CN,Ut)) + 2 A: + iil^ (t) P{h) . (5.5) 

We still need to determine (CAf,ki(f) CAf,k2(^); which requires that we use quantum 
perturbation theory. 

Interestingly, the contribution from the field-redefinition is the same as the non- 
Gaussianity predicted by the "5A^" formalism [Q, 27, |2^, this connection was originally 
established by Seery and Lidsey |30|, who noted that that the field redefinition used in the 
bispectrum calculation is similar to the transformation from comoving to fiat gauge, and 



^(^ always goes to zero at late times since ^ oc /a^^Q far after horizon crossing [13| 
^°FT[/(x)](k) is defined as 



FT[/(x)](k)^ / dV(x)e-*-^ 
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Arroja and Koyama |31, 32], who observed that the gauge transformation on large scales 
is given by the "5A^" formalism. To verify this result, note that in the "(5A^" formalism, 
one supposes that (/> is completely Gaussian and that any non-Gaussianity is generated by 
non-linearities in the relationship between and </>: 



to second order in 84), this gives Q = Ql — ( ^^d4>'^ ) ^i' "^^^^^ Cl = —1/2 {d\uH / d(j)')~^ 6cj} 
is linear in 5(f) and Gaussian. Then, 

{C') = -{2nf 6('\ Ek.) 2 [Pdki)Pdk2) + Pdk2)Pdh) + Pdk3)Pdh)] , (5-7) 

which gives us, for single-field inflation, 

(C^>.3«.... = -(2-)^^(^)(Ek.) Pc(^i)^c(^3) 2 + ; (5.8) 

this is the same as the bispectrum contribution from the second term in (^]^) after we 
correlate with (^ks- 

5.3 Quantizing (^^ 

Here, we quantize Cn, following the methods of [^]. We start by expanding (j\f: 

C7v(x) = I a.ke^'^-'' • (5.9) 



From (5J.), we find the equation of motion for Cn,\s.' 



dt 

To quantize Cat, we promote ^Af.k to an operator 

C,N,\<i{t) = Uk{t)ak + ul{t)ai^^ . (5.11) 

Ok) Ok annihilation, creation operators for ^]\f, respectively, which obey the canonical 
commutation relation [ak,o|^,] = (27r)^(5^^)(k — k'). Uk,ul. are two independent solutions 
to the equation of motion ( ^.lOD for (^^r, normalized to reflect the appropriate vacuum 
condition for large k; Uk,ul. are known as the "mode functions" of ^Af- Note that since 
is the same as ^5 to lowest order in (, the mode functions for ^tv are the same as those for 
(s and for 

At this point, we can find the power spectrum of Cn (which is the same as the power 
spectrum of ("): 

(CkCk'> = (2^)='|tXfc|V3)(k + k'), 

so 

P{k) = \uk\^. 
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5.4 Applying the in-in formalism 

The in-in formalism p^ , 35 sometimes called the Keldysh-Schwinger formalism, has 



become an important tool in the calculation of cosmic correlation functions |14, 30, 37, 
38|, |39|, |3l|, It specifies the expectation value of quantum operators at some time 



given a known input state in the past. Here, we are interested in calculating (Civ(i)CAf(O) 
given that the universe was in a vacuum state in the far past. 

The question of selecting the proper initial vacuum state is non-trivial and has been 
written about extensively |33, 42, ^]^^. The eff'ect of vacuum state choice is usually 
accounted for by mixing the mode functions u, u* to form new mode functions u,u*. In 
this section, we do not explicitly specify mode functions, so we are not forced to choose 
an initial state. In the following section, where we work examples, we choose the standard 
Bunch-Davies vacuum, which assumes that the mode functions of fields limit toward their 
standard Minkowski-space form at early times and/or deep inside the horizon. 

The formula for the tree level in-in formalism (sufficient for our calculation) is 

{C\t)) = -i f dt'(0|[C2(t),ifKO]|0>, (5.12) 

where Hj is the interaction Hamiltonian in the interaction picture and (o| . . . |o) means 
that we take the expectation value in the free-field theory. In principle, there is a distinc- 
tion between variables in the interaction picture and in the Heisenberg picture but this 



can largely be ignored without any complications (see |43] for more about the interaction 



picture). A different, but relevant, subtlety is that the time integral contour in ( 5.12 ) picks 
up a small imaginary part at early times, so that we have written the lower limit of the 
integral as —(1 — ie)oo. This contour should be familiar to those who have done stan- 
dard perturbation calculations in QFT; in both cases, the purpose of the imaginary part 
is to perform the calculation in the perturbed (rather than the free) vacuum. If the mode 
functions are in the Bunch-Davies vacuum, this usually makes the integrand exponentially 
decay at early times so that the integral converges. 

It is usually true, including in this case, that Hj = —Lint. Then, inserting our expres- 
sion (5^) for Lint into the formula above (see appendix^ for the explicit calculation), we 
arrive at: 



where 



K = iUf,^{ri) I dr] 

-oo(l — ie) 



1 "^0 2 /*2/ \ , ^ "^0 2)2 \ , 

■a u^, (r?) + ^j^a kiUk^{r]) + 



2H^ 



3 d / (/>0 1 '/'o I 



H'^ dt \ 6nH 



+ c.c. . 



(5.13) 



35 1 has a list of other references on the in-in formahsm. 

37 1 has a hst of other references on the effect and intricacies of initial vacuum selection. 
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Using this previous formula in ( |5.5D and correlating with Cks) we get: 

(C=^> = (2vr)3 5('\ Ek.) P{h) |p(fci) 2 + + . (5.14) 

This is the main result of the paper: a formula for the bispectrum in the squeezed limit, 
given in terms of the mode functions u,u* . Note again that no slow-roll approximation has 
been used. 



6. Checking the consistency relation in specific cases 
6.1 Slow-roll inflation 

Slow-roll inflation is the simplest model of inflation and, additionally, is entirely compat- 
ible with current observational data. It corresponds to a nearly de Sitter (dS) expansion 
produced by a single inflaton field. We assume H = — (1/2)(/)q is small and so is (pQ. More 
precisely, we define dimensionless slow-roll parameters 

to be small. We then perform all our calculations to lowest order in these parameters. 

In dS space, we can always choose % and such that a = —(r]H)^^. This is still true 
for slow-roll inflation, at least to lowest order in slow-roll. 

The equation of motion ( 5.1C| ) for and thus for u becomes: 



dt 
so that 



d (a^ ilk) 2 

— ak Uk = 0, 



0Q V 



and 



,2 ^f* 1 



Spectral tilt It is reasonable to suppose that in P{k), H'^/cpQ = H'^it^) / (f)'^{t^), where 
t=K is the time that the mode k crosses the horizon, i.e. k w a{t^)H. Then, dk/dt^ w 
a{t^)H^ ^ kH, so d{lnk) = Hdt^. Then, 

4HPmn ^ 1 dlngy^g ^ jk_ ,6.3) 

In A; H dt^ H(t)Q ' 
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Bispectrum in the squeezed limit Using ( ^.14 ): 

^ ^ V 'J-oo(l-ie) ^ ' 



--P{ki) 



H'^ ki 



/)g 2A:i j_oo(i-ie) 



-1\2, 



+ C.C. 



Then, since we take the — t- hmit, to determine the bispectrum at the present we let 



2iJ4 V i 



i,2 2 V2iA;i 



if2 1 



1 1 

kiT] 2i 



2kir] 



+ c.c. 



ri=ri 



(6.4) 



So, 



(C^^) = (2vr)3 5(3)( ^k.) P{k^)P{k,) 2{^ + %) 

= (27r)3 5^% YMi) P{ki)P{h) (1 - ns), (6.5) 
matching the result of |14] and verifying the consistency relation. 
6.2 Power-law inflation 

As another example, consider power-law inflation (where the consistency relation has not 
previously been verified). It provides an opportunity to test the consistency relation in a 
fully non slow-roll situation. 

In power-law inflation [U, |4|, we suppose F(0o) = ge~^'^° , where g and A are 
dimensionless parameters. We can exactly solve the scalar field equations ( [4.2[ ) to find 



In 



3-e 



and 



H 



A2 



we have defined the dimensionless parameter 

H _ 

~ 16^' 

which has some similarities to the slow-roll parameter esR except that here, we only require 
e < 1. 
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Then, we can choose ao and to such that 

a = t^'\ 

It is useful to define the quantity 

3 e 



2 ' 1-e' 



then, slow-roll inflation corresponds to ~ 3/2. 
We can solve (5.10) for u 



Ukiv) 



c e 



e V 1 - e 



1 

" l-e 



exp 



where we have defined c = -\^/A{2nf/'^ and Hi^\z) is a Hankel function; the relevant 
properties of Hankel functions are listed in appendix 
Thus, we find that 



P{k) = \u{k)W^^^ 



c ^T{v)f ( e \--(2 
1 e2 7r2 \l-e \k 



2v 



(6.6) 



and 



\-ns 



din 



din /c 



2i/-3 



2e 
l-e' 



(6.7) 



We can also calculate that {l/2)4>'^/ H"^ = e and 4>o/(I)qH = — e. Then, ( 5.14 ) becomes 

^ J-oo(l-ie) 



Observe that 



c e 



e V 1 - e 



c / e 



exp 



vH^^\-kr^) + {-kri) 



dHi^\-kT]) 
d {—kr]) 



e V 1 - e 



exp(!|^ + ^]fe(-,r//W(-M, 



so that 



K = -2c^P{ki) { i 



oo(l— ie) 



dz z 



Hi'\{z) > F(2)(. 



+ c.c. 



where z = —iki-q. 

Using the equations from appendix we find 

f"Oo(l— ie) 



dz z 



Hi'\z) 



+ c.c. = u 
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so 

„2 



K = -2c'P{ki) { --{v - 1) + --v\ = P{h] 



vr vr J 1 — e 

Then, as desired, 

2e 



= {27rf ^k,) P{k,)P{k,) (1 - ns). (6.9) 
6.3 An exactly scale-invariant model by Starobinsky 

It's worth noting that, in the previous two models, the slow-roU parameters e and fj were 
constant, which greatly simplified our calculation of K because we could pull out the factors 
of e in the first two terms of ( ^.13 ) and because the last term is exactly zero^^. Motivated 



to find an infiation model that predicts a scale-invariant spectrum, Starobinsky in 2005 
described another exactly solvable inflation model with varying (not necessarily small) 
slow-roll parameters [^]; this model has 1 — = 0, so that, by the consistency relation, 
it predicts an exactly vanishing squeezed-limit bispectrum and is therefore an interesting 
test case. To characterize the model, define z = acp/H (for this section, we will drop the 
subscript on </>o). His model is essentially specified by assuming z to be of the form 

B 

V ' 

where 5 is a positive free-parameter. This is enough to show that^^, 

1 e''"'! (\ 



''' = BWk[~k^'' 

For canonical single- field models, the criterion for a mode to freeze outside the horizon is 
that k"^ <C z" /z because then = 1/z^ = H"^ /a?(j)^ For this model, this implies that 
/c^ ^ 2/r/^, so that the natural choice for the end of inflation is when all modes are frozen 
outside the horizon, i.e. that rj = 0.^^ 

If we suppose that inflation ends when r/ = 0, P{k) = l/2B^k^ which is exactly 
scale-invariant, as desired. 

For this model, 

V{(t>) = 3H{4>f - 2H'i^), 

where 

H{cl)) = ^e^'l^ e-^'/U4> + C j , (6.10) 



^^Because 4>o/<}}oH + (l/2)(^g/H^ = -v + 2e. 

^^Note that our definition of u appears different from ust in [|l^ because 1) u — ust/z, and 2) we shifted 
the phase of u so that -P(fci) ~ ^^(end of inflation). 

^^In principle, one might object that e rises above unity slightly before rj reaches 0, implying that inflation 
ends earlier. However, it's not clear that e is a better indicator of whether the horizon size is shrinking. 
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and C is another free parameter. For C = 0, there is an infinite period of inflation (i.e. 
d, d > 0) that becomes slow-roll as i;^ — >• oo. For C small and negative, there is a finite 
period of inflation which gets arbitrarily long as C — )• 0. We will thus focus on these cases. 

For this model, it is sometimes useful to use (j) as the dynamical variable instead of ij; 
for the cases of interest, there is a one-to-one correspondence between these variables. We 
will indicate derivatives with respect to (p using hats (e.g. H = dH/d(f)). We flnd, for C 
small and negative. 



1 



B 
2ao 



-ao- 



74 



BH 



2ao- 



2e-02/4 



where ao = a[(j) = 0). For reference, note that 



+ 



2 4>H 4^2 



-1 



2 H 

bW 



e-^^/^d^ + C 



(hi) 277^ 
7] rj"^ 



.11 



(6.12) 



For this model, K becomes a fairly involved integral, which is best performed in terms 
of (j). The integral in K can be evaluated after performing a (somewhat complicated) series 
of integrations by part: 



drj 



2^4« %i +2^4« ^l«/c,+2^2« at 



2H^ 




g2jfci77 
,2ikiTi 



rff] (4 



2kir] 



V 



2f 

^2 



There is a slight subtlety when evaluating this quantity at rj — t- — oo(1 — ie), for the lower 
bound of integration in K, because in the Starobinsky model, rj{(f)) doesn't limit to — oo, 
and in fact, the model behaves unusually for (p somewhat less than zero. However, one can 
imagine that for ~ 0, the model merges with another inflation model or a static universe. 
Then, we just let rj go to —oo and the e^'^^^ factor will make the integrand go to zero. 
Thus, 



-oo(l— ie) 

and the term in brackets in ( |5.14| ) is 



drj ... + c.c. 



1 



B'^^l 



1 



(pfl 
V 



2fj 



a2 



\uk,?2 



(pH 



-rf + 'qcpi) + 2f/^) 



P{ki 



B^k^ 



/2 
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where cf)cnd = 't'i'n = 0) </'end = 2erfc ^— -^j, erfc is the complementary error 

function), yielding 

/nl = {B^kj/al)e-'f''nJ\ (6.13) 

Since one expects to find /nl = 1 — = for the Starobinsky model, this result 
initially appears confusing. However, as we will explain, /nl is non-zero only to the extent 
that some of the conditions of the consistency relation are violated. First, we will discuss 
the factor B'^kf/ag. As mentioned earlier, the comoving horizon size is approximately 
1/(1/77) = ?? so that, at = 0, the horizon size is \r]{(p = 0)| w B/2ao, since C is taken to 
be small. Thus, the ratio Bki/aQ specifies how close a mode ki is to the horizon size at the 
start of inflation. Since (for C ^ 0) the Starobinsky model has a finite inflationary period, 
this ratio also determines how long the mode will spend outside the horizon; note that, 
formally, all modes cross the horizon as — )• but some of them will spend only extremely 
short times outside the horizon. The second factor e~'^end/2 has a similarly straightforward 
interpretation: as C — t- 0, (j)end 00 and the inflationary period gets arbitrarily long, so all 
modes spend an increased time outside the horizon. Both of these factors are reasonable 
since the consistency relation should only hold if the perturbations are outside the horizon 
"long enough" to freeze out and to stop varying; our result is interesting in that it shows 
quantitatively how /nl approaches the consistency relation prediction as the parameters 
approach the conditions it assumes. 



7. Discussion 



In this paper, we developed a technique for calculating squeezed-limit correlation functions. 
We applied it to the single-field bispectrum with a canonical kinetic term and arrived at 
( 5.14 ), a formula for the squeezed-limit bispectrum which relies on no slow-roll approxi- 
mation. We first verified that our technique matched the known squeezed-limit result for 
slow-roll inflation. We then performed the calculation for power-law inflation to arbitrary 
order in slow-roll and explicitly verified that it matched the consistency relation prediction. 
Finally, in the Starobinsky model, we again saw that the consistency relation held, though 
we found corrections that depended on the ratio of the size of the ki mode and the horizon 
size and on the length of the inflationary period. Thus, one has to be careful that the 
assumptions of the consistency relation are satisfied before expecting it to hold. Even for 
a mode that leaves the horizon, if the subsequent inflationary period is sufficiently short, 
it could have an enhanced squeezed-limit bispectrum. 

Our technique can also be applied to cases with non-standard kinetic terms. This was 
done in [p^ j; following our method, the author derived the formula for the squeezed-limit 
bispectrum for any single- field model where the action for the scalar is given by P{X,(f)), 
where X = —{l/2)g^^'^d^4'dy(j). He verified the validity of the consistency relation for 
power-law inflation where the speed of sound 7^ 1, as well as for second-order slow-roll 
inflation. 

Besides providing explicit verification for the consistency relation, our technique pro- 
vides a method of calculating (CkiCk2)^j^ i an important part of the consistency relation's 
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derivation. Furthermore, as compared to the proof of the consistency relation, it is done 
purely in Fourier space, i.e. it provides a complementary approach. 

We would also like to point out that our formalism could be useful for calculating other 
squeeze-limit correlation functions, since our primary assumption is that wavemodes which 
stretch far enough outside the horizon turn classical. In particular, our technique could be 
adapted for calculating the squeezed limit trispectrum or multi-field inflation bispectrum. 
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A. Some explicit calculations 



A.l Field redefinition calculation from subsection 5.2 



If Cs(t,x) = Civ(t,x) +/3(t)CL(t,x)Cjv(t,x), then 



(C5,ki(t)C5,k2(0> = 

= I d^xid^X2 (Cs(xi)Cs(x2))e-^(''-''i+'^^-''^)« 

~ (C7V,kiC^,k2> / d^Xid^X2 (C7V(X1)C7V(X2)) [/5 Cl(xi) + /3 Cl(x2) 



-j(ki-xi+k2-X2) 



(CAf,kiCiV,k2) + 

+ 2 /5 / n (^^'^l'^'^2 a,q3 (CjV.q,lCjV.q,2)^ € 
' =(27r)3FT(gi)5(qi+q2) 

(C7V,kiC^,k2) +2/3(t)CLM+k2(t) FT [{e{Ax,t))]{h) 



«((q3+qi-ki)-xi + (q2-k2)-X2) 



matching the result from earlier. As noted earlier in a footnote, FT[/(x)](k) is defined as 



FT[/(x)](k)^ J dV(x) 



-ikx 
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A. 2 In-in calculation from subsection 

First, let us calculate from (5.4), using Q and 



Hi{t) = -L 



int 



a,- 



qi-q2 - 



(2^) 

1 <^r3(qi-q2)' 



H rrffl" 

16i/6 



Cqi Cq2 + tt9 



+ 



Cqi Cq2 



(A.l) 



where we have neglected the C,NdiC,Ndid '^C.l term because it will not contribute in the 
squeezed limit (this is further discussed at the end of this subsection). 
Also, 

(0|[Cki(t)Ck2(t),Ck,(t)Ck2(i)]|0> = 




<\2 



i2TTf (u{ki,t)u{k2,t){u(qi,t')uiq2,t'))* - c.c.) x 
X [6{ki + qi)(5(k2 + q2) + S{ki + q2)5(k2 + qi)] 



Uk, (t) flki + (t) a^ki ) (""te (t) Oka + <2 (t) 0^2 

X (^Ug, (t) + (t) alq^ ) (^^52 (t) + U*^ (t) al 

(^Ug^ (t ) + n*^ (t ) alq^ ) (^Ugj (t ) + (t ) alq^ ) x 



(A.2) 



There are also variants where the (" terms with q indices have time derivatives, for example 
where Cqi(*) is replaced by Cqi(0; ™ these cases, we simply apply the time derivative to 
the corresponding mode function, so that tiqi(t) nqi(t). 
Then, 



(CAf.kiCAf.ka) 



^ ,.d^qid^q2. 

to 

X 



(27r)e 



-qi-q2 • 



4;^(i')«(i')qi • q2 (0 |[Ck,Ck2,CqiCq2]l 0) + 

1 



Cki Ck2 ) Cqi Cq2 



+ 
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1 



^0 



1 i4 



(0(o 



Cki Ck2 ) Cqi Cq2 



(t')a(t')ki • k2 u{ki,t)u{k2,t){u{ki,t')u{k2,t'))* 
.,^(kl•k2) 



1 



+ 



H 2m 



X 2n(A;i , t)u{k2,t) {u{ki,t')u{k2,t')y 



c.c. 



kiRika 



CL,ki+k2-f^ 



(A.3) 



where K is defined in ( 5.14| ). 

Now let us return to the piece we neglected earlier: 



Lintit) C 



^0 . 



2H^ 



a%s^iCs^^^-\L. 



(A.4) 



This contributes 



Hi{t) C 



^3 



(27r)6 2i?4 



q2 • (-qi - q2) 



qi-q2 



-qi - q2) 



CqiCq 



2 SqiSq2 ' 



(A.5) 



so that 



(CkiCka) C 



(pQ „3 q2 • i-qi - q2. , ^ 



Cki Ck2 ) Cqi C 



q2 



to 



2H^ 



(kl+k2)2 



u{ki,t)u{k2,t) [u{ki,t')u{k2,t')y 



c.c. 



+ ki o k2 



j'dt'^{t')a{t'fCLM+^2 [u\ki,t) {u{k,,t')u{ki,t')) 



c.c. 



(A.6) 



However, this is the same contribution we would get from a term 



4iJ4 



(A.7) 



which is clearly non-dominant in the squeezed limit. 
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B. Properties of Hankel functions used in subsection |6.2 



Here, we collect the properties of Hankel functions that are used in our calculations (these 
can be found in any standard reference, e.g. |^]). [Equations that do not specify H^^^ or 
H^"^^ apply to both and arguments shown as z can be complex while arguments shown as 
X are real.] 

HW(x) = Hi'>{x), (B.l) 
z^H,{z) - vH,{z) = zH,+i{z), (B.2) 

J dxx[H,{x)f = y [[H,{z)f - H,^i{x)H,+i{x)\ (B.3) 
iz) ^ ^ (f ) ^ ± (^) ^ real x (B.4) 

H^^){z) —?■ J —e^''^^^^^^^) large \z\, |argz| < tt (B.5) 

V TTZ 

References 

[ll A. A. Starobinsky, Spectrum of relict gravitational radiation and the early state of the 
universe, JETP Lett. 30 (1979) 682-685. 

K. Sato, First Order Phase Transition of a Vacuum and Expansion of the Universe, Mon. 
Not. Roy. Astron. Soc. 195 (1981) 467-479. 

A. H. Guth, The Inflationary Universe: A Possible Solution to the Horizon and Flatness 
Problems, Phys. Rev. D23 (1981) 347-356. 

A. D. Linde, A New Inflationary Universe Scenario: A Possible Solution of the Horizon, 
Flatness, Homogeneity, Isotropy and Primordial Monopole Problems, Phys. Lett. B108 
(1982) 389-393. 

A. Albrecht and P. J. Steinhardt, Cosmology for Grand Unified Theories with Radiatively 
Induced Symmetry Breaking, Phys. Rev. Lett. 48 (1982) 1220-1223. 

V. F. Mukhanov and G. V. Chibisov, Quantum Fluctuation and Nonsingular Universe. (In 
Russian), JETP Lett. 33 (1981) 532-535. 

S. W. Hawking, The Development of Irregularities in a Single Bubble Inflationary Universe, 
Phys. Lett. B115 (1982) 295. 

A. A. Starobinsky, Dynamics of Phase Transition in the New Inflationary Universe Scenario 
and Generation of Perturbations, Phys. Lett. B117 (1982) 175-178. 

A. H. Guth and S. Y. Pi, Fluctuations in the New Inflationary Universe, Phys. Rev. Lett. 49 
(1982) 1110-1113. 

flOl J. M. Bardeen, P. J. Steinhardt, and M. S. Turner, Spontaneous Greation of Almost Scale - 
Free Density Perturbations in an Inflationary Universe, Phys. Rev. D28 (1983) 679. 

[Ill W. Fischlcr, B. Ratra, and L. Susskind, Quantum Mechanics of Inflation, Nucl. Phys. B259 
(1985) 730. 



- 21 - 



H. Kodama and M. Sasaki, Cosmological Perturbation Theory, Prog. Theor. Phys. Suppl. 78 
(1984) 1-166. 

V. F. Mukhanov, H. A. Feldman, and R. H. Brandenberger, Theory of cosmological 
perturbations. Part 1. Classical perturbations. Part 2. Quantum theory o J perturbations. Part 
3. Extensions, Phys. Kept. 215 (1992) 203-333. 

J. M. Maldacena, Non-Gaussian features of primordial fluctuations in single field inflationary 



models, JHEP 05 (2003) 013, |astro-ph/0210603 |. 

E. Komatsu and D. N. Spergel, Acoustic signatures in the primary microwave background 



bispectrum, Phys. Rev. D63 (2001) 063002, |astro-ph/0005036 |. 

E. Komatsu et. at, Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) 



Observations: Cosmological Interpretation, arXiv : 1001 . 4538. 

P. Creminelli and M. Zaldarriaga, Single field consistency relation for the 3-point function, 



JCAP 0410 (2004) 006, iastro-ph/0407059 



A. A. Starobinsky, Infiaton field potential producing the exactly flat spectrum of adiabatic 



perturbations, JETP Lett. 82 (2005) 169-173, |astro-ph/0507193 



C. Kiefer, D. Polarski, and A. A. Starobinsky, Quantum-to-classical transition for 



fluctuations in the early universe. Int. J. Mod. Phys. D7 (1998) 455-462, [ [gr-qc/9802003| ] 
E. Komatsu, The Pursuit of N on- Gaussian Fluctuations in the Cosmic Microwave 



Background. PhD thesis, Tohoku University, 2002. astro-ph/0206039 



A. A. Starobinsky, Stochastic De Sitter (Inflationary) Stage in the Early Universe, in Field 
Theory, Quantum Gravity and Strings (H. De Vega and N. Sanchez, eds.), April, 1986. 
107-126. 

S. Habib, Stochastic inflation: The Quantum phase space approach, Phys. Rev. D46 (1992) 



2408-2427, |lgr-qc/9208006 



S. Rcnaux-Petel, On the squeezed limit of the bispectrum in general single field inflation, 
arXiv : 1008 ■ 026C| . 



R. L. Arnowitt, S. Deser, and C. W. Misner, The dynamics of general relativity, in 
Gravitation: an introduction to current research (L. Witten, ed.), ch. 7, pp. 227-265. Wiley, 
New York, 1962. 

D. Wands, K. A. Mahk, D. H. Lyth, and A. R. Liddle, A new approach to the evolution of 
cosmological perturbations on large scales, Phys. Rev. D62 (2000) 043527, 



| astro-ph/0003278 |. 

D. H. Lyth and Y. Rodriguez, Non-gaussianity from the second-order cosmological 



perturbation, Phys. Rev. D71 (2005) 123508, [|astro-ph/0502578|] 



A. A. Starobinsky, Multicomponent de Sitter (Inflationary) Stages and the Generation of 
Perturbations, JETP Lett. 42 (1985) 152-155. 

D. S. Salopek and J. R. Bond, Nonlinear evolution of long wavelength metric fluctuations in 
inflationary models, Phys. Rev. D42 (1990) 3936-3962. 

M. Sasaki and E. D. Stewart, A General analytic formula for the spectral index of the density 
perturbations produced during inflation. Prog. Theor. Phys. 95 (1996) 71-78, 



astro-ph/9507001 



- 22 - 



D. Seery and J. E. Lidsey, Primordial non-gaussianities from multiple-field inflation, JCAP 



0509 (2005) Oil, iastro-ph/0506056 |. 

F. Arroja and K. Koyama, Non-gaussianity from the trispectrum in general single field 



inflation, Phys. Rev. D77 (2008) 083517, ]arXiv : 0802. 1167 1. 

K. Koyama, Non-Gaussianity of quantum fields during inflation, Class. Quant. Grav. 27 



(2010) 124001, I arXiv: 1002.06001 



N. D. Birrell and P. C. W. Davies, Quantum Fields in Gurved Space (Cambridge Monographs 
on Mathematical Physics). Cambridge University Press, April, 1984. 

J. S. Schwinger, Brownian motion of a quantum oscillator, J. Math. Phys. 2 (1961) 407-432. 
S. Weinberg, Quantum contributions to cosmological correlations, Phys. Rev. D72 (2005) 



043514, |hep-th/0506236| 



D. Seery and J. E. Lidsey, Primordial non-gaussianities in single field inflation, JCAP 0506 
(2005) 003, | |astro-ph/050369^ . 



X. Chen, M.-x. Huang, S. Kachru, and G. Shin, Observational signatures and 
non-Caussianities of general single field inflation, JCAP 0701 (2007) 002, | hep-th/ 06050451 . 

D. Seery, J. E. Lidsey, and M. S. Sloth, The inflationary trispectrum, JCAP 0701 (2007) 027, 



| astro-ph/0610210 |. 

D. Seery, M. S. Sloth, and F. Vernizzi, Inflationary trispectrum from graviton exchange, 
.JCAP 0903 (2009) 018, 



arXiv: 0811. 3934 



X. Chen, B. Hu, M.-x. Huang, G. Shiu, and Y. Wang, Large Primordial Trispectra in General 



Single Field Inflation, JCAP 0908 (2009) 008, JarXi v: 0905. 3494 1. 

F. Arroja, S. Mizuno, K. Koyama, and T. Tanaka, On the full trispectrum in single field 



DBI-inflation, Phys. Rev. D80 (2009) 043527, |arXiv: 0905. 3641 



B. Greene, M. Parikh, and J. P. van der Schaar, Universal correction to the inflationary 
vacuum, JHEP 04 (2006) 057, [ |iep-tli/0512243| ]. 

S. Weinberg, The Quantum Theory of Fields, Volume I: Foundations. Cambridge University 
Press, May, 1995. 

L. F. Abbott and M. B. Wise, Constraints on Ceneralized Inflationary Cosmologies, Nucl. 
Phys. B244 (1984) 541-548. 

F. Lucchin and S. Matarrese, Kinematical Properties of Ceneralized Inflation, Phys. Lett. 
B164 (1985) 282. 

D. H. Lyth and E. D. Stewart, Constraining the inflationary energy scale from axion 
cosmology, Phys. Lett. B283 (1992) 189-193. 

S. Weinberg, Cosmology. Oxfort Univ. Pr., Oxford UK, 2008. 

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products. Academic Press, 
7 ed., March, 2007. (A. Jeffrey and D. ZwiUinger, eds.). 



-23- 



